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In this paper a new theory for the calculation^ the radiation resistance of antennas in gyroelectric 
media is presented. -This new theory-ensures the irreversibility and finiteness of the radiation 
resistance. 



1. Introduction 

In recent years many investigators [e.g., Staras, 1964; Weil and Walsh, 1965] have contributed 
significantly to the problem of calculating the radiation resistance of a dipole antenna immersed 
in an electrically anisotropic homogeneous medium whose anisotropy is due to a magnetostatic 
biasing field B . However, the method they have used to calculate the radiation resistance — the 
so-called conventional method — leads in general to an infinite value for the radiation resistance, 
and on this account the method and the result it yields are open to question. Some have inter- 
preted this infinity as a type of resonance in the medium, while others have attributed it to the 
infinitesimal size of the source. It is our contention that the difficulty is of a-more basic nature and 
is due not simply to the size of the source but to the method of calculation. In this paper we show 
that when the radiation resistance is calculated in accord with the laws of reversibility and irre- 
versibility, the radiation resistance turns out to have a physically reasonable value. Since radia- 
tion resistance is on the same footing as circuit resistance in the sense that they both are measures 
of irreversible power, we require that only the irreversible part of the power be used in calculating 
radiation resistance. 

Our point of departure is the conventional expression for the time-average power P. Using 
combinations of outgoing and incoming fields we split P into its irreversible part P irr and its re- 
versible part P rev such that P = Pi rr + P re v and show that P irr is free from singularity whereas P rev 
is not. Thus, for the radiation resistance R we obtain R = CP lTT where C is a constant. This result 
obeys the required conditions of finiteness and irreversibility. Previous investigators have used 
the expression R = CP = C{P irr -\- P rey ) and therein lies the source of their difficulty because this 
expression, due to the presence of P rev , violates not only the required irreversibility but also the 
condition of finiteness. We present a general formula for the irreversible power radiated by an 
arbitrary current distribution varying harmonically in time. We show that this formula reduces 
to the conventional result when B = (isotropic media) or when B ( > = ^ (uniaxial media). 

To construct an expression for the irreversible part of the power emitted by a source we recall 
that in the case of an accelerating point electron in vacuum the combination of half the retarded 
minus half the advanced field is free from singularity, [Dirac, 1938], and corresponds to the irrever- 
sible power radiated by the electron [Schwinger, 1949]. We extend this idea of taking a combina- 
tion field to the case of a monochromatic source radiating into a lossless, homogeneous, gyroelec- 
tric medium. 
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2. Formulation of the Problem 

The radiation resistance R of an antenna is defined by the relation 

R = CP ivr (1) 

where C is the reciprocal of the squared amplitude II* of the current Re Ie~ l0)t feeding the antenna 
and P irr is the irreversible part of the time-average real power P emitted by the antenna. Accord- 
ingly, the problem of calculating radiation resistance amounts to the problem of calculating from a 
knowledge of the antenna current density Re J(r)e~ la)t the irreversible part of P. 
It is well known that the power P is given by 



| Re j J*(r) • E 0Ut (r)dV (2) 



2 

where the integration extends throughout the volume occupied by the antenna current and E out 
denotes the electric vector of the outgoing wave generated by the antenna. However, it has not 
been recognized that P generally has not only an irreversible part P irr but also a reversible part 
P rev as displayed by the relation 

r / irr ' * rev {*>) 

Indeed, the question of reversibility and irreversibility has never before been raised in this connec- 
tion. This is perhaps due to the fact that the dual nature of P may be overlooked in certain special 
cases without affecting the value of the radiation resistance. For example, when the antenna is 
in vacuum or in a simple lossless unbounded medium P rev is identically zero and hence the conven- 
tional definition R = CP is indistinguishable from our definition (1). On the other hand, in the 
case where the ambient medium is a lossless, homogeneous, unbounded gyroelectric medium — the 
very case we are presently interested in — P rev is not necessarily identically zero and consequently 
the operation of finding the irreversible part of P is a crucial step in the calculation. 

Our problem, therefore, is to extract from P an explicit expression for P irr involving the an- 
tenna current density J(r) and the dyadic Green's function of the gyroelectric medium. 

3. Time-Reversal Transformation 

Since the anisotropy of the gyroelectric medium which surrounds the antenna is produced by 
a magnetostatic biasing field, the reversal of time must be accompanied by the reversal of the 
biasing field if Maxwell's equations are to be covariant under the time-reversal transformation. 
Formally, the time-reversal transformation may be thought of as a transformation from a ref- 
erence frame K with space coordinates x, y, z, time coordinate t, and biasing field B , to a reference 
frame K' with the same space coordinates, but with time coordinate t' ', and biasing field B^. In 
transforming from K to K' the space coordinates remain unchanged whereas the time and the 
biasing field transform according to the rule 

t^t'=-t, Bo^B; = -Bo- (4) 

Since the microscopic Maxwefl-Lorentz equations are known to be covariant under this trans- 
formation, we demand that the macroscopic Maxwell's equations for a lossless gyroelectric medium 
be likewise covariant. Using the convention that primed quantities are referred to K' and un- 
primed ones to K, we display this covariance by stating that in K we have 

V X H(r) - J(r) - io)e(Bo) ■ E(r), V • e(B ) ■ E(r) = p(r) 

V X E(r) = ia)/ioH(r), V • H(r) = (5) 
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and in K' we have 

VxH'(r) = J'(r)-/'co6(B ( ;) • E'(r), V- e(Bi) • E'(r) = p'(r) 

VxE'(r) = ^ H'(r), V- H'(r) = 0. (6) 

It follows from the covariance of these macroscopic Maxwell equations that the field and source 
quantities in K are related to those in K' by the transformations: 

E(r) = E'*(r), H(r)=-H'*(r), J( r ) = -J'*( r ) 

p(r) = p'*(r), e(B ) = e*(B ). (7) 

We can verify that the above transformation law for the dielectric tensor is indeed obeyed when 
Bq = — Bo, by recalling that the dielectric tensor of a lossless gyroelectric medium has the property 

e(B ) = e*(-B ). 

The physical significance of these transformations (7) may be more easily appreciated in the 
time domain. For example, in the case of the electric vector we have 

E'(r, t')=Re E'(r)e-^'=Re£*(r) e i(tit = Re E(r) e-^ f = E(r, *) (8) 

or 

E'(r,-*) = E(r,t) (9) 

which shows that the past and the future of an electric field are the same. Similarly the past 
and the future of any of the source or field quantities (except for minus signs) are the same. It 
is also revealing to consider Poynting's vector. If by S' and S we denote the complex Poynting 
vector in K' and K, respectively, we have the transformation rule 

S' = \ E'(r) x H'*(r) = - |e» x H(r) =- S* (10) 

whose real part yields 

ReS' = -ReS. (11) 

Thus we see that the real part of Poynting's vector changes sign under time-reversal. 

Up to the present point of our discussion it appears that in a lossless gyroelectric medium the 
electromagnetic field and the real part of the complex Poynting's vector are completely reversible. 
To answer the question of where the irreversibility comes from we must remember that the fields 
we have been examining are total fields and that Maxwell's equations actually permit two in- 
dependent wave solutions, one representing an outgoing wave E out , and the other an incoming wave 
E in . In the time domain the two independent solutions are the retarded and the advanced wave, 
which are related to E out and E in by 

E rot (r,0 = ReE out (r)e-^ 

E adv (r,0 = ReE in (r)e--< (12) 
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in K, and by 

E; t (M') = ReE; ut (r)e-'-' 

E; dv (r,«') = ReE;(r)e-"' ,13) 

in K' . Since E ret (r, t)= Ea dv (r, £'), it follows from (12) and (13) and similar equations for the 
magnetic fields that 

E 0Ut (r)=E; n *(r) 

E to (r) = E *(r) 

H out (r) = -H i ' n *(r) 

H ln (r) = -H *(r). (14) 

Relations (14) show that an outgoing wave in K transforms into an incoming wave in K' and vice 
versa. 

To interpret the result in (14) in terms of Poynting's vector we denote the complex Poynting 
vectors of the outgoing and incoming waves in K by S out and S in , and in K' by S ( ' )Ut and Sj' n , respec- 
tively, and note that the transformations 

Sout = \ E„ut(r) X H* ut (r) =- | E,'*(r) X H; n = - S,' n * 

S in = |E^r)XH* n (r) = - |E' * ut (r) X H' ut (r) = -S' * t (15) 

lead to 

Re(S 0Ut + S in ) = Re(- Si* ± S;*) = ± Re(S ( ; ut T S/ n ). (16) 

This equation demonstrates that the combination Re(S 0Ut — S in ) is irreversible and Re(S 0Ut 4- S in ) 
is reversible under the time reversal transformation (4). 

From the complex Poynting's vector theorem, which in the present case of a lossless gyro- 
electric medium has the form 

V • (E X H*)=- J* • E + ico^oH • H* - io)E • e(B ) • E*, (17) 

we get 

Re V • (E X H*) =- Re J* • E (18) 

which on integration yields 



Re [n-Sda = - | Re f J*(r) • E(r) dV 



(19) 



where the volume integration extends throughout the volume occupied by the current and the 
surface integral extends over any closed surface with outward normal n, enclosing the current. 

1316 



Applying this relation (19) to S out and S in we obtain the relations 

Re I n • S out da = - - Re I J*(r) • E out (r) dV 

Re f n • S in da = - | Re J J*(r) * E in (r) dF (20) 

whose difference and sum yield 

Re J n • (S out + S in ) da=- | Re | J*(r) • [E out (r) + E in (r)] dF. (21) 

Since by (16) Re(S 0Ut — S in ) is irreversible and Re(S 0Ut -h S in ) is reversible, we see from (21) that 

Re f J*(r) • [E out (r) - E in (r)] dV (22) 

J I 

expresses irreversible power and 



/ 



Re J J*(r) • [E out (r) + E in (r)] dV (23) 

expresses reversible power. The irreversibility of (22) and the reversibility of (23) can also be 
made evident by submitting them to the field transformations of (14) and the current density trans- 
formation of (7). 

Recalling expression (2) for P and using the algebraic identity 

E ou t = 2^ (,ut_ E in ) + -(E ()Ut + E in ) (24) 

we find that 

P=-|Re[j*(r)-[E ut(r)-E ln (r)]dF 

-| Re [ J*(r) • [E out (r) + E ln (r)] dV. (25) 

Thus we split P into its irreversible part (first term on right) and its reversible part (second term 
on right). 

Consequently we can write 

* = "irr "I" 'rev (26) 

where P irr and P rev are given by 

P m =-\ Re J J*(r) • [E out (r) - E ln (r)] dV (27) 

P rev =-\Rej J*(r) • [E 0Ut (r) + E in (r)]dF (28) 

or equivalently by 

Pirr = | Re | n • (S out - S in ) da (29) 

^rev = | Re J n • (Sout + Sib) da. (30) 

To find P irv we must use either expression (27) or expression (29). These expressions are note- 
worthy because they place in evidence the necessity of including the incoming field or the incoming 
Poynting's vector. 
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4. Regular Behavior of the Difference Field 

Since we have 

P = P irr + P rev (31) 

and since Pis not free from singularity [Staras, 1964; Weil and Walsh, 1965], the question of how 
the possible singularity in P is distributed between the terms P irr and P rev arises. Clearly, to 
answer this query one examines expressions (27) and (28) for P irr and P rey . Here we make the 
reasonable assumption that the volume integral of J itself is bounded, and therefore the deter- 
mining quantities to be examined are the difference field E out — E in , and the sum field E out + E in . 
We know that E out and E in satisfy the same differential equation, i.e., 

V X V X E 0Ut -(D 2 )LC() e • E out = i(ofi J (32) 

Vx VxE in -^ ( ,€-E in = ^«,J. (33) 

However, at infinity the boundary condition on E out is different from the boundary condition on 
E in , and therefore E out and E in cannot be the same function. 

Taking the difference of (32) and (33) we obtain the homogeneous equation for the difference 
field 

VxVx(E out -E in )-G>Vo€-(E out -E in ) = (34) 

which holds everywhere. This equation implies that E out — E in is free from singularity. Using 
this result and recalling that E out ¥^ E in , we see from (27) that P irr is finite. 

On the other hand, taking the sum of (32) and (33) we see that the sum field satisfies the in- 
homogeneous equation 

V x V x (E out + E in )-a> 2 A*o e • (E out + E in ) = 2 iw/xoJ (35) 

which exhibits the singular behavior of the sum field. Since it is the sum field that carries the 
singularity, we see from (28) that the singularity appears in P rev - 



5. Irreversible Power in Terms of the Dyadic Green's Function 

To reduce expression (27) for the irreversible power to a simpler computational form we intro- 
duce the dyadic Green's functions which relate E out and E in to the source current density J. 
From the linearity of Maxwell's equations we can write 



(36) 

(37) 



E out (r, ± Bo) = i(ofio J T out(r, r', ± B„) ■ J(r') dV 

E in (r,±Bo)=iwAo J r ln (r,r\±Bo)-J(r')dK' 
where T out and T in are the dyadic Green's functions which satisfy 

V X V X f out (± Bo)- wVo6(± Bo) • T out (± B ) = u8(r- r') (38) 

VX Vxr in (±Bo)-co 2 jUoe(±B„)-r in (±B )=u6(r-r') (39) 
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where u=unit dyadic. Although T out and T in satisfy the same differential equation, they are no! 
identically equal because they must satisfy different boundary conditions at infinity. Since 
€(— B()) = €*(Bo) we find by taking the conjugate complex of (39) thai 

V X V X r* F Bo) - wVo€(± B () ) • r*(+ Bo) = uo(r - r '). (40) 

Comparing (38) with (40) we see that T out (±B ) and T ^(+B () ) satisfy the same equation. In the 
special case where B o = we know that T out = Tf n , and thus we conclude that 

r* n (+Bo) = r ()Ut (±Bo) (41) 

holds for all values of B . 

With the aid of relation (41) we obtain from (36) and (37) the following expression for the 
difference field: 

E out (r, Bo) - E in (r, B ) = ico^ ) I [ r out (i\ r', B„) - r * ut (r, r\- B«)] • J(r')dV . (42) 

Substituting this expression into (27) we see that the irreversible power is given by 

Pirr= £ f 5 Im[j*(r) • [r„u<<r.r\ B„)r *„(r, r'.-B,,)] ■ }(r')dVdV. (43) 

To simplify further this relation we note that 

r ut(Bo)-r* ut (-Bo) = Re[r ou ,(Bo)-r ou t(-Bo)]+i [m[r out (B )+r out (-B )] (44) 

and accordingly write P lrr in the following form 

P, rr =^Re[j*(r)-Im[r ou ,(r,r',Bo) + r out (r,r',-Bo)]-J(r')drar 

+ ^ lm|j*(r)-Re[r out (r,r',Bo)-r out (r.r'.-Bo)]-J(r')r/Fjr. (45) 

This is the desired expression for the irreversible power. It shows that from a knowledge of J 
and r out we can find / J irr , and from a knowledge of P irr we can find in turn the radiation resistance 
by the formula R = CP xrT . 

For the sake of completeness, we also include the expression for the reversible power: 

P rev = ^Im J J*(r) • Re [r^ 

+ ^Re| J*(r) ' Im [IWr, r', B )- r out (r, r',-B )] • J(r') dVdV . (46) 

6. Dipole Antenna 

Let us now consider the case of a dipole antenna in a gyroelectric medium. For a dipole 
antenna situated at the origin of the coordinates we have 

J(r) = -icop8(r). (47) 
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Substituting this expression into (45) we see that the irreversible power radiated by the dipole is 
given by 

Pirr = ^p-Im[r ut(0,0,Bo)+r out (0,0,-Bo)]-p (48) 

which is clearly an even function of the biasing field B as it should be. 

If we had used expression (2), as one would do in the conventional method, we would have 
obtained 

P=^ p .Imr out (0,0,Bo)-p. (49) 

Evidently expressions (48) and (49) are different. However, in the extreme cases where B = 
(isotropic media) or where Bo = °o (uniaxial crystals), the conventional expression (49) becomes 
identically equal to our expression (48). Hence, when B = oo or when B () = 0, it is permissible 
to use expression (49) for the calculation of the radiation resistance [Kogelnik, 1960; Kuehl, 1962], 
but otherwise it is not. When B has a finite value, expression (49) cannot be used for this purpose 
because it does not represent irreversible power. When B is finite one must use expression (48). 

7. Conclusions 

From the above discussion we can draw the following conclusions. For an antenna in a gyro- 
electric medium, the time average power P emitted by the antenna has an irreversible part P irr 
and a reversible part P rex . The irreversible part of the power does not change sign under the time- 
reversal transformation and hence represents the time-average real power that is absorbed by the 
sphere at infinity; the reversible part of the power, on the other hand, changes sign under the time- 
reversal transformation and accordingly has the nature of reactive power. When the biasing field 
Bo is zero or infinite P rev is identically zero, and there is no difference between P and P irr . How- 
ever, for finite values of B () there is a profound difference. This means that for finite values of 
Bo the radiation resistance R, which is a measure of the power absorbed by the sphere at infinity, 
must involve only the irreversible part of the power. 

Although the resulting electromagnetic field of an antenna in a gyroelectric medium consists 
of only the outwardly moving wave, it is expedient to introduce the inwardly moving wave because 
the difference of the outgoing and incoming waves yields P irr and their sum yields P rev . More- 
over, consideration of the difference and sum fields shows that P irr is free from singularity, whereas 
P rev is not. 

By means of the dyadic Green's function the explicit representations of P irr and F rev take a 
particularly simple form. Indeed, they require only a knowledge of two dyadic Green's functions, 
one of which can be obtained from the other by reversing the sign of the biasing field. 
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